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Abstract 

We present a simple nonlinear relaxation equation which contains 
the Debye equation as a particular case. The suggested relaxation 
equation results in power-law decay of fluctuations. This equation 
contains a parameter defining the frequency dependence of the di- 
electric permittivity similarly to the well-known one-parameter phe- 
nomenological equations of Cole-Cole, Davidson-Cole and Kohlrausch- 
Williams- Watts. Unlike these models, the obtained dielectric permit- 
tivity (i) obeys to the Kramers-Kronig relation; (ii) has proper be- 
haviour at large frequency; (iii) its imaginary part, conductivity, shows 
a power-law frequency dependence cr ~ cj" where n < 1 corresponds 
to empirical Jonscher's universal relaxation law while n > 1 is also ob- 
served in several experiments. The nonlinear equation proposed may 
be useful in various fields of relaxation theory. 

The frequency dependence of the dielectric permittivity s{u) plays a cru- 
cial role in numerous aspects of condensed matter science. Among the great 
number of publications, we can mention the textbooks [H [2] and review arti- 
cles [31 H] . Classical theory of polar dielectric relaxation, which is connected 
with the e{co), is based on the hnear Debye equation 

Here W is the angular distribution of dipoles, Wq is its equilibrium value 
and r is a relaxation time. The solution of the Eq. ([1]) describes an expo- 
nential relaxation, which results in the Debye dielectric permittivity 

= £00 + ez)(w), eD(w) = 7^— (2) 
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where £00 is the dielectric permittivity at a; ^ 1/r and Eg = £{0). However, 
the dependence occurs quite seldom in experiments [H El El S] • Expla- 
nation of a non-Debye dispersion e{uj) is carried out in frames of the models 
which can be tentatively divided into three directions, (i) Models which 
assume the relaxation time r to be a stochastic variable. Taking different 
distributions for r one can obtain from the equation different e{uj) depen- 
dencies, including those matching the experimentally observed curves (see, 
for instance, Refs. [HI [6]). Apparently this theoretical approach to explain 
the experimental data is not faultless, (ii) Electrical network model used 
for the microstructural media [7j. (iii) Models which consider influence of 
manybody effects on the relaxation [U SI El E] ■ The model suggested here is 
close to the latter type. 

The kinetic equation for spatial distribution of dipoles includes a colli- 
sion integral which model is right hand side of Eq. ([1]) (see, for instance, Ref. 
|10j). In general kinetic theory this integral is a nonlinear functional of the 
distribution. It is therefore natural to take the right hand side of the relax- 
ation equation as a nonlinear function. The simplest nonlinear generalization 
of the Debye equation ([T|) is the following relaxation equation: 

dW _ {W - Wo)^+i 
dt ~ T ■ 

The parameter q in the equation defines the relaxation law, shape of 
absorption line and dispersion dependence for the real part of e. In this regard 
the parameter q is analogous to the parameters of Cole-Cole, Davidson-Cole, 
Kohlrausch- Williams- Watts phenomenological curves |£j . In the present case 
the fitting parameter is introduced at an earlier stage, i. e., into a differ- 
ential relaxation equation. This leads to a number of principal differences 
comparing to the aforementioned phenomenological curves. 

A nonlinear relaxation equation 

dw _ - W^^^ 

dt T ' 

which is similar to the Eq. (E|), also includes the only shape parameter q, 
though it leads to more sophisticated mathematical formulae. 

Obviously, for noninteger q the right hand side of the Eq. (jSj) displays 
nonanalytic (in the sense of Cauchy-Reimann conditions) behaviour which 
cannot follow from kinetic theory. For small deviations from the equilibrium, 

U = W -Wo, (4) 

the right hand side of the relaxation equation should have Debye form ([1]). 
However, this is a rather difficult question whether U can be considered as 
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Figure 1: Hypothetical view of the colhsion integral I{U) = U + U"^ (solid 
curve) and approximation of it in the Debye linear model I{U) ~ U (dashed 
curve) and in the proposed nonlinear model I{U) ^ 2U^'^ (points). Multiplier 
2 in the latter equation defines the difference between the relaxation times 
in the linear and nonlinear models. 



small enough. We assume that the U values relevant to the experiment are 
not small so that the relaxation equation can be approximated by the simple 
formula (see Fig. [T]). 

From the Eq. one easily obtains the relaxation function 

f/(^) = (g^/r + C)-^/^ (5) 

where q is positive number to satisfy the condition W — > Wq at t — > oo. The 
constant C = Uq'^ ensures the initial condition U{0) = Uq. Then 

f/(t) = f/o(l + gtf/o7r)-^/^. (6) 

Function coincides with the Debye relaxation function U {t) = Uq exp(— t/r) 
for g ^ + 0. 

The general relation between the electric field E{t) and the electric dis- 
placement field D{t) will be used. The principle of causality gives for the 
isotropic media [11] 

D{t) = e^E{t) + f Kit - e)E{e)de, (7) 
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where /t may differ from tlie relaxation function U{t) only by a constant 
positive multiplier 



«:(t) = A(l + ^/r*)-l^ (8) 
and an effective relaxation time 

is introduced. 

Multiplier A can be calculated for q < 1 from the equation ([7j) with the 
dc field E. Taking into account the equation D = SgE for this case, 

^= ^^(£s-£oo), q<l. (10) 
qr* 

For g > 1 the field E cannot be assumed as a dc since integral in the 
equation ([7]) diverges. Multiplier A cannot be calculated in this case. This 
problem is discussed below. 

To determine an ac permittivity £{oo) let us consider a periodic field 
E{t) = Eoe-"^^ in the Eq. d?]). Then 

e(u) = eoo + e(u), (11) 



e = 6^+162 = A (1 + t/r*)-^^'^e'^'dt 
Jo 



= — (-za;r*)i/V~*""*r (1 - 1/q, -iut*) , 

CO 

where r(-, ■) is an incomplete F-f unction \13\. In a more common form one 
can write e{uj) as 

7 A 

e = ti-it2 = («u;r*)^/«e*""*F (1 - iujt*) . (12) 

It is easy to see the dependence e on tu is determined by effective relaxation 
time r*. The actual relaxation time, r, which enters the Eqs. ([5]), ([6]), can be 
measured by the methods of time domain spectroscopy [HI [15]. Since e(c<;) 
is Fourier-Laplace transform of a smooth function it obeys to the Kramers- 
Kronig relation [12]. Function e{uj) has an integrable singularity at w = 
when q > 1 and this singularity is an inaccuracy of the model. In fact, it 
is well known that behaviour of a function in frequency domain near w = 
is determined by the function in time domain at t — )• oo. However, for very 
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large t the function W is close to the function Wq and equation becomes 
incorrect. So the equation f lT^ cannot be used with g > 1 in the narrow 
frequency region near u = 0. 

Imaginary part of the permittivity, €2, is connected with the conductivity 
of a matter, a, by the well-known equation 

= ^^2M (13) 

and ciuj) has no singularity at a; = for any q. Its behaviour at small cu is 
the following 

r £T{1 - l/g) cos(7r/2g) (ur*)'/", q > 1/2 

a(w) ~ <^ A(a;T-)2 ^ 1 /o (14) 



A(a;T'') 
47r(l-2g) 



, g < 1/2 



It is interesting to note that the dependence cr ~ a;^, which follows from 
the Eq. (fT^ for q < 1/2, coincides with its dependence in Debye theory 
when g = 0. The dependence a ^ u"^ with n < 1 has been observed in 
many experiments (Jonscher's universal relaxation law [21 El [ISl [IZ] ) . Such 
a dependence follows from the equation f|T^ for g > 1. Dependence a ~ 
with n > 1 which follows from the equation f|T^ for 1/2 < g < 1 has 
been also observed in some types of materials, e. g., in glassy 0.3(xLi2O.(l- 
x)Li2O)0.7B2O3 [IB], in mixed compounds of (NH4)3H(S04) 1.42 (Se04)o.58 [IS] 
and K3H(Se04)2 single crystals [SO]- Sometimes at very low frequencies, 
u < uq, in both (n < 1 and n > 1) cases the frequency dependence of 
conductivity is weak, plateau-like. It does not contradict the proposed model 
since in the low-frequency region the conductivity may be connected with 
other processes, magnitude of which becomes equal to ([H]) when a; ~ cjq [S] . 

For — )■ 00 

that coincides with the asymptotic behavior of e/j ([2]). We remind that the 
dependence ei ~ o;"^ is a general feature of a dielectric permittivity at large 
frequencies [TT]. Nevertheless, for many materials the decrease of ei with 
the increasing frequency is slower than u~'^. It may be connected with the 
influence of faster relaxation processes which induce extra high-frequency 
maxima on the dielectric losses curve. 

Clearly, the suggested model requires a detailed experimental verification. 
Here we demonstrate that it may be adjusted with the well-known empirical 
dependencies, e. g., with the Cole-Cole equation 

Ae ^ , , 

''''' = TTJ^^^ A. = ..-.^. (16) 
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Figure 2: Imaginary part of the dielectric permittivity, €2, calculated accord- 
ing to the equation (|T2|) (solid curve) and to the Cole-Cole equation (fT6|) 
(dashed curve). All parameters are given in the text. 

Figure [2] shows the frequency dependence of €2 for a = 1/2. The parameters 
of nonlinear relaxation are chosen as g = 0.732 and r* = A.GGtcc- Such 
a large difference of r* and tcc implies a shift of dielectric losses curve 
maximum to the low-frequency region. 

In conclusion, the nonlinear relaxation equation is proposed. The equa- 
tion leads to a power decay of the fluctuations. Frequency dependencies of 
the dielectric permittivity and dielectric losses are in agreement with the well- 
known empirical dependencies. Derived relaxation equation may be useful 
to describe different relaxation processes. If a relaxation of the dimensional 
values is under consideration, U in should be equal to W/Wq — 1. 
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